Search for: All records

Abstract Subduction zones are host to some of the largest and most devastating geohazards on Earth. The magnitude of these hazards is often measured by the amount of energy they release over short periods of time, which itself depends on how much stored energy is available for the geologic processes that drive these hazards. By considering the energy transfer among processes within subduction zones, we can identify the energy inputs and outputs to the system and estimate the stored energy. Due to the multiscale nature of subduction zone processes, developing an energy budget of subduction zone hazards requires integrating a wide range of geologic and geophysical field, laboratory, and modeling studies. We present a framework for developing mechanical energy budgets of upper crustal deformation that considers processes within the magmatic system, at the subduction zone interface, distributed and localized deformation between the arc and trench, and surface processes that erode, transport, and store sediments. The subduction energy budget framework provides a way to integrate data and model results to explore interactions between diverse processes. Because fault mechanics, sediment transport and magmatic processes within subduction zones do not act in isolation, we gain insights by considering the common energetic elements of the subduction zone system. Building energy budgets reveals gaps in our understanding of subduction zone processes, and thus highlights opportunities for new interdisciplinary research on subduction zone processes that can inform hazard potential. 

40Ar/39Ar detrital sanidine (DS) dating of river terraces provides new insights into the evolution and bedrock incision history of the San Juan River, a major tributary of the Colorado River, USA, at the million-year time scale. We dated terrace flights from the San Juan−Colorado River confluence to the San Juan Rocky Mountains. We report >5700 40Ar/ 39Ar dates on single DS grains from axial river facies within several meters above the straths of 30 individual terraces; these yielded ∼2.5% young (<2 Ma) grains that constrain maximum depositional ages (MDAs) and minimum incision rates. The most common young grains were from known caldera eruptions: 0.63 Ma grains derived from the Yellowstone Lava Creek B eruption, and 1.23 Ma and 1.62 Ma grains derived from two Jemez Mountains eruptions in New Mexico. Agreement of a DS-derived MDA age with a refined cosmogenic burial age from Bluff, Utah, indicates that the DS MDA closely approximates the true depositional age in some cases. In a given reach, terraces with ca. 0.6 Ma grains are commonly about half as high above the river as those with ca. 1.2 Ma grains, suggesting that the formation of the terrace flights likely tracks near-steady bedrock incision over the past 1.2 Ma. Longitudinal profile analysis of the San Juan River system shows variation in area-normalized along-stream gradients: a steeper (ksn = 150) reach near the confluence with the Colorado River, a shallower gradient (ksn = 70) in the central Colorado Plateau, and steeper (ksn = 150) channels in the upper Animas River basin. These reaches all show steady bedrock incision, but rates vary by >100 m/Ma, with 247 m/Ma at the San Juan−Colorado River confluence, 120−164 m/Ma across the core of the Colorado Plateau, and 263 m/Ma in the upper Animas River area of the San Juan Mountains. The combined dataset suggests that the San Juan River system is actively adjusting to base-level fall at the Colorado River confluence and to the uplift of the San Juan Mountains headwaters relative to the core of the Colorado Plateau. These fluvial adjustments are attributed to ongoing mantle-driven differential epeirogenic uplift that is shaping the San Juan River system as well as rivers and landscapes elsewhere in the western United States. 

Volcanic provinces are among the most active but least well understood landscapes on Earth. Here, we show that the central Cascade arc, USA, exhibits systematic spatial covariation of topography and hydrology that are linked to aging volcanic bedrock, suggesting systematic controls on landscape evolution. At the Cascade crest, a locus of Quaternary volcanism, water circulates deeply through the upper $\sim$1 km of crust but transitions to shallow and dominantly horizontal flow as rocks age away from the arc front. We argue that this spatial pattern reflects a temporal state shift in the deep Critical Zone. Chemical weathering at depth, surface particulate deposition, and tectonic forcing drive landscapes away from an initial state with minimal topographic dissection, large vertical hydraulic conductivity, abundant lakes, and muted hydrographs toward a state of deep fluvial dissection, small vertical hydraulic conductivity, few lakes, and flashy hydrographs. This state shift has major implications for regional water resources. Drill hole temperature profiles imply at least $81$km ${}^3$of active groundwater currently stored at the Cascade Range crest, with discharge variability a strong function of bedrock age. Deeply circulating groundwater also impacts volcanism, and Holocene High Cascades eruptions reflect explosive magma–water interactions that increase regional volcanic hazard potential. We propose that a Critical Zone state shift drives volcanic landscape evolution in wet climates and represents a framework for understanding interconnected solid earth dynamics and climate in these terrains. 

Linear feature analysis plays a fundamental role in geospatial applications, from detecting infrastructure networks to characterizing geological formations. In this paper, we introduce linkinglines, an open-source Python package tailored for the clustering and feature extrac- tion of linear structures in geospatial data. Our package leverages the Hough Transform, commonly used in image processing, performs clustering of line segments in the Hough Space, and then provides unique feature extraction methods and visualization. linkinglines em- powers researchers, data scientists, and analysts across diverse domains to efficiently process, understand, and extract valuable insights from linear features, contributing to more informed decision-making and enhanced data-driven exploration. We have used linkinglines to map dike swarms with thousands of segments associated with Large Igneous Provinces in Kubo Hutchison et al. (2023). 

Abstract Large igneous provinces erupt highly reactive, predominantly basaltic lavas onto Earth’s surface, which should boost the weathering flux leading to long-term CO₂drawdown and cooling following cessation of volcanism. However, throughout Earth’s geological history, the aftermaths of multiple Phanerozoic large igneous provinces are marked by unexpectedly protracted climatic warming and delayed biotic recovery lasting millions of years beyond the most voluminous phases of extrusive volcanism. Here we conduct geodynamic modelling of mantle melting and thermomechanical modelling of magma transport to show that rheologic feedbacks in the crust can throttle eruption rates despite continued melt generation and CO₂supply. Our results demonstrate how the mantle-derived flux of CO₂to the atmosphere during large igneous provinces can decouple from rates of surface volcanism, representing an important flux driving long-term climate. Climate–biogeochemical modelling spanning intervals with temporally calibrated palaeoclimate data further shows how accounting for this non-eruptive cryptic CO₂can help reconcile the life cycle of large igneous provinces with climate disruption and recovery during the Permian–Triassic, Mid-Miocene and other critical moments in Earth’s climate history. These findings underscore the key role that outgassing from intrusive magmas plays in modulating our planet’s surface environment. 

Article Published: 27 May 2024 Explosive 2018 eruptions at Kīlauea driven by a collapse-induced stomp-rocket mechanism Josh Crozier, Josef Dufek, Leif Karlstrom, Kyle R. Anderson, Ryan Cahalan, Weston Thelen, Mary Benage & Chao Liang Nature Geoscience volume 17, pages572–578 (2024)Cite this article 1357 Accesses 430 Altmetric Metricsdetails Abstract Explosive volcanic eruptions produce hazardous atmospheric plumes composed of tephra particles, hot gas and entrained air. Such eruptions are generally driven by magmatic fragmentation or steam expansion. However, an eruption mechanism outside this phreatic–magmatic spectrum was suggested by a sequence of 12 explosive eruptions in May 2018 at Kīlauea, Hawaii, that occurred during the early stages of caldera collapse and produced atmospheric plumes reaching 8 km above the vent. Here we use seismic inversions for reservoir pressure as a source condition for three-dimensional simulations of transient multiphase eruptive plume ascent through a conduit and stratified atmosphere. We compare the simulations with conduit ascent times inferred from seismic and infrasound data, and with plume heights from radar data. We find that the plumes are consistent with eruptions caused by a stomp-rocket mechanism involving the abrupt subsidence of reservoir roof rock that increased pressure in the underlying magma reservoir. In our model, the reservoir was overlain by a pocket of accumulated high-temperature magmatic gas and lithic debris, which were driven through a conduit approximately 600 m long to erupt particles at rates of around 3,000 m3 s−1. Our results reveal a distinct collapse-driven type of eruption and provide a framework for integrating diverse geophysical and atmospheric data with simulations to gain a better understanding of unsteady explosive eruptions. 

Along subduction zones, high-relief topography is associated with sustained volcanism parallel to the plate margin. However, the relationship between magmatism and mountain building in arcs is poorly understood. Here, we study patterns of surface deformation and correlated fluvial knickpoints in the Columbia River Gorge to link long-term magmatism to the uplift and ensuing topographic development of the Cascade Range. An upwarped paleochannel exposed in the walls of the Gorge constrains unsteady deep magma flux, the ratio of intrusive to extrusive magmatic contributions to topography, and the impact of magmatism on Co- lumbia River incision since 3.5 million years ago. Geophysical data indicate that deep magma influx beneath the arc axis is ongoing and not aligned with the current locations of volcanic edifices, representing a broad regional influence on arc construction. 

Combining visual and sonic representations of data can make science more accessible and help reveal subtle details. The recent decade-long eruption of Hawaii’s Kīlauea Volcano offers a prime example. 

We develop a method for estimation of parameters of an elastic plate resting on a Winkler-type elastic foundation solely from data on the vertical displacements of the plate. The method allows one to estimate components of the external body force density field, plate thickness, elastic foundation stiffness parameters, horizontal displacements of the plate, and stresses. The key idea of the method is that multiple plate models are used simultaneously, namely the proposed reduced three-dimensional (R3D) plate model, the Mindlin plate model, and the thin plate model. The three plate models form a hierarchy of elastic plate models based on assumptions imposed on stresses, with the R3D plate model being the most generalized model and the thin plate model being the most constrained one. The hierarchical relationship among the plate models allows one to incorporate prior information into the estimation technique. The applicability of the proposed estimation method is illustrated by a numerical example.